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In this talk I review measurements and calculations of the flux of neutrinos pro- 
duced by interactions of cosmic rays in the atmosphere. The main reason for 
interest in this subject is the apparent anomaly between the predicted and the 
observed ratio of Ve to Vf^. With the advent of Super-Kamiokande, we are on the 
threshold of an order of magnitude increase in the amount of data available to 
study the problem. My goal in this talk is to describe the current status of the 
subject, both for contained events and for neutrino-induced upward muons. 



1 Introduction 

Because of their small cross sections neutrinos were the last component of the 
secondary cosmic radiation to be measured, although they are the most nu- 
merous particles in the GeV energy range at sea level. Markov El suggested 
how upward and horizontal muons deep underground could be used as a sig- 
nal of high energy neutrinos, and GreisenEI described a neutrino detector like 
the modjsrn water detectors in which neutrino interactions could be observed 
directlyu Neutrino-induced horizontal muons at the predicted level j^ere ob- 
served a few years later in deep mines in IndiaQ and in South Africau 

The deep detectors built to search for proton decay have now accumu- 
lated more thaii, a thousand contained events. The large water Cherenkov 
detectors, 1MB □ and KamiokandeP dominate the statistics. Although the 
total events rates are consistent with the expectation for interactions of atmo- 
spheric neutrinos, the ratio of electron-type to muon-type neutriiLQs is signif- 
icantly higher than predicted for the water Cherenkov detectorsQ^Q The first 
hint of this anomaly came already in 1986 with the_pbservation by 1MB of 
fewer than expected muon decays among their eventsB The most well-defined 
class of events is the contained single-ring events. These are mostly charged- 
current quasi-elastic events (e.g. -\- n ^ p -\- fj,~) with with an admixture of 
neutral-current events in which a single pion is produced. A recent statement 
of the anomaly for < GeV neutrinos at Kamiokande isEl 
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for contained single-ring events. Tlae,icnn tracking calorimeters, hmsseiier, find 
results consistent with no anomaly or with a smaller anomalyHlS 
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The rate of interactions for neutrinos inside a detector of mass M (in g) 

Rate = NaM J dE^ j dEi j dQ ME,,,fl) ^ e{Ee), (1) 

where A^^ is Avogadro's number. The factors in the integrand are the dif- 
ferential flux of Vi; the differential cross section to produce the corresponding 
lepton, i; and the efficiency, e, for its identification and detection. The expres- 
sion for the rate of neutrino-induced muons is similar except that the target 
mass is 

M{n) - R^{E^) X Aifl), (2) 

where R^{E^) is the muon range and A{^) is the projected area of the detector 
as seen from the direction Q. = {6, </>). 

Explanations for the anomaly in the contained events have been sought 
in all three factors of Eq. |^. I will discuss calculation of the flux of at- 
mospheric neutrinos separately in the next section. Several approximations 
have been made in the treatment of the cross sectionSpEpr neutrino interac- 
tions in oxygen inside the water detectors. Engel et alEJ conclude, however, 
that the neglected physics cannot account for the anomalous /i-to-e ratio ob- 
served at Kamiokande and 1MB. The reason is that the lepton momenta are 



high enough that corrections a: 
same way. Beam tests at KEK 



muons and electrons in very nearly the 
confirm the efficiencies determined from 
simulations for misidentifying muons as electrons and vice versa. 

Another possibility is that there is some contamination of events that are 
not due to interactions of atmospheric-neutrinos. The suggestion that there 
was an excess of events due to p — > i/DecIi isjaresumably eliminated by the fact 
that the anomaly persists at higher energyES RyazhskayaEj has suggested that 
extra electron-like events are really cascades from neutral pions produced inside 
the detector by interactions of entering neutrons generated outside the detector 
in interactions aLjatmospheric muons in the rock. Kamiokande argues against 
this by showingliB that they have a tt*^ production rate consistent with neutral 
current production by neutrinos at the level expected and already accounted 
for in their simulations Soudan argues against this explanation by analysis of 
their shield eventsEjj 



Such possibi 
Various analyses 



s highlight the importance of internal checks on the data. 
I agree that the Kamiokande and 1MB data are consistent 
with each other after accounting for differences in exnosure, geomagnetic cut- 
off and energy threshold, although Beier and Frank Ej note a possible small 
discrepancy in the electron spectra at low momentum. There is a hint of the 
expected smaller interaction rate at the 1MB detector during the period of 
maximum solar activity between 1989 and 1991 as compared to the earlier 
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part of their data collection from 1986-88, a period of minimum solar activity. 
One expects fewer events during solar maximum when the low energy primary 
cosmic rays are partially excluded from the inner solar system. The effect 
should be more noticeable in the overall rate at 1MB than at Kamiokande be- 
cause of the higher local geomagnetic cutoff at Kamiokande, which excludes 
a large fraction of the lower energy primaries in any case. During the first 
period, (exposure of 3.4 kT yrs) 1MB found 236 events (139 e-like and 97 fj,- 
like)t3 The corresponding numbers for the last 4.3 kT yrs of the 7.7 kT-yr. 
exposure wereQ 271 (186 e-like and 85 /i-like). The total expected by simple 
extrapolation of the first period would be 297 ± 20. The difference (297 — 271), 
though not statistically significant, is about what would be expected due to 
solar modulation effects. The fact that the decrease shows up only in the //-like 
events is not consistertLwith solar modulation. Presumably it is an accident of 
low statistics. StanevcJ has emphasized the importance of using the expected 
solar cycle variations as a probe of atmospheric neutrino data, which will be 
possible with larger data samples. 
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Before turning to the discussion of the neutrino fluxes, I display the in- 
tegrand of Eq. ^ Fig. 1 shows the "response-curves" for observed neutrino- 
induced muons. What is pfotted is the rate of muons integrated over muon 
energy per logarithmic interval of neutrino energy. The four classes^of events 
are fully-contained interactions, Kamiokande "multi-GeV" eventsp£l entering 
muons thai-Stop in the detector and neutrino-induced throughgoing muons. 
The figure l3 is specific to the Kamiokande detector, but similar distributions 
could be constructed for any detector. The curves rise from low energy as 
the neutrino cross section increases with energy. For i/-induced muons the ef- 
fective volume also increases with energy as the range increases. Eventually 
the growth of range and cross section slow and the steep cosmic-ray spectrum 
(together with pion interaction) cuts off the signal of atmospheric neutrinos at 
high energy. The response curves are useful when considering possible expla- 
nations of the flavor anomaly in terms of neutrino oscillations. 



2 Flux of atmospheric neutrinos 

The analysis of atmospheric neutrino expeiimeiits has depended mainly on four 
calculations of atmospheric_neutrinosE30cj'L3 The jmlculations of G. Barr, 
Gaisser aniStanev (BGS)@, Honda et al. (HKHM)E3 and Bugaev and Nau- 
mov (BN) cB are completely independent of each other. The neutrino flavor 
ratio is the same within 5% in all these calculations, but there are some signif- 
icant differences in normalization and shape of the calculated neutrino energy 
spectra. Many sources of uncertainty cancel in the calculation of the ratio of 
Ve/vfi- Thus the theoretical uncertainties are much smaller in the ratio than in 
the normalization. Fogli and Lisi lHI have shown how to make the comparison 
between expected and measured neutrino interactions in this situation. 

Suzuki i has compared the measured spectra of electrons and muons in 
single ring neutrino interactions with full simulations starting from the fluxes 
of HKHM, BGS and BN. The calculation of BGSEj gives the steepest spectrum, 
predicting more muons with momenta below 600 MeV/c than observed, but 
agreeing withJJjie measured spectrum of electrons. In contrast, the neutrino 
spectra of BNc3 are nearly in agreement with the muon spectrum down to 2QQ 
MeV/c but predict fewer electrons than observed—The calculation of HKHMtZl 
is intermediate but closer to the results of BGSEj. 

Table 10 compares the neutrino fluxes of the three calculations. The first 
part of the table shows the neutrino spectra separately for Ve and i'^ in three 
energy intervals, normalized to BGS= 1.00. The second part of the table shows 
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Table 1 : Comparison of calculated neutrinos fluxes at Kamioka 





0.4-1 1-2 2-3 


Ve. + Ve 

0.4-1 1-2 2-3 


BGS 

HKHM 

BN 


1.00 1.00 1.00 
0.90 0.95 1.04 
0.63 0.79 0.95 


1.00 1.00 1.00 
0.87 0.91 0.97 
0.62 0.74 0.87 





0.4 < < 1 GeV 


BGS 

HKHM 

BN 


0.99 0.89 0.49 
0.99 0.84 0.48 
0.98 0.76 0.50 



the neutrino ratios in the energy interval between 0.4 and 1 GeV. Here 

^e/^i — 1 



to reflect the smaller interaction cross section for antineutrinos. This crucial 
ratio isjaearly the same in all cases. 

Wecil have investigated the sources of difference among the calculations by 
substituting one-by-one different sets of assumptions from^tlie various papers 
into the framework of the BGS calculation. In that workcS the spectrum of 
neutrinos Vi was expressed as a convolution of the primary cosmic-ray spec- 
trum, the geomagnetic cutoffs and the yield per nucleon of Vi. 

4>^^{n) = (t>p® Rp®Yp^^^ (3) 

In this equation Rp[Ra) is the geomagnetic cutoff for protons (nuclei) incident 
on the atmosphere from the direction — {0^(j)\. The three terms on the 
right-hand-side represent respectively neutrinos from primary hydrogen, from 
protons bound in incident nuclei and from neutrons bound in incident nuclei. 
The separation is necessary because neutrino production depends on energy- 
per-nucleon of the incident cosmic rays but the geomagnetic cutoff depends on 
magnetic rigidity (gyroradius) . Nuclei and protons of the same energy per nit 
cleon differ by a factor of A/Z in magnetic rigidity. To a good approximation's 
the yields of neutrinos from nuclei can be calculated as if the incident nucleons 
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were unbound. (This approximation somewhat overestimates the production 
of neutrinos from pions produced in the target fragmentation regiDu for the 
fraction of the flux due to incident nuclei. It is not used in Ref. Ej.) In the 
energy region important for contained events, approximately 80% of the neu- 
trinos are produced by cosmic-ray hydrogen (free protons) and most of the rest 
come from helium nuclei. 

The form of the BGS calculation (Eq. |) makes it possible to trace the 
effects of different assumptions through the calculation. As an example, con- 
sider the 1 GeV neutrino flux at Kamiokande at solar minimum. (Kamioka has 
the highest cutoff for downward cosmic rays of the nucleon decay detectors, so 
the effect of differences in cutoff is maximum. Comparison at solar minimum 
maximizes the effects of diffeiences in assumed prirnary spectrum.) We com- 
pare the calculation of BGS@ with that of HKHMB. The treatment of the 
geomagnetic cutoff in BGS neglected the "penumbra" effect. A more accurate 
treatment of the cutoffs in by HKHM reduces the 1 GeV flux by a factor 0.87. 
The primary spectrum assumed in HKHM is higher than in BGS, which gives 
a factor 1.25. Yields of GeV neutrinos are approximately 15% lower in HKHM, 
giving a factor of 0.85. The product 0.87 x 1.25 x 0.85 = 0.93 gives a net 7% 
lower GeV neutrino flux in HKHM than in BGS. 

The main result of our comparison^ islkat Jie biggest source of difference 
among the three independent calculations ESEZIo is the treatment of produc- 
tion of low energy pions in collisions of 10 to 30 GeV protons with nuclei of the 
atmosphere. In the BN calculation the inclusive cross sections for production 
of < 3 GeV pions is significantly lower than in the calculations of BGS and 
HKHM. In fact, it is quite similar to the spectrum of pions in proton-proton 
interactions. In contrast, both HKHM and BGS use representations of pion 
production that give significantly more low energy pions in collisions on nitro- 
gen than for pp collisions. It is this feature of the BN calculation that gives 
rise to their characteristically harder neutrino spectra with relatively few low 
energy (< 2 GeV) neutrinos. At higher energy the neutrino fluxes of the three 
calculations are in better agreement. 

Measurements of rnupns at high altitude can be used to constrain the neu- 
trino spectra. Perkins lj has calculated neutrino spectra starting from mea- 
sured muon spectra, including preliminary results of the MASStj experiment. 
He concludes that the higher flux is preferred. Several groups have now mea- 
sured the muon flux during the ascent of their Jmlloon-borne detectors. The 
spectra of negative muons first reported in Ref . CJ have now been publishedE3 
The IMAX-experiment has given a preliminary report of their measurements 
of muons E3, sjoA the HEAT experiment is in the process of analyzing their 
measurements t3 
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Both the Japanese groupEj and we have now pubhshed calculations of 
the neutrino fluxea^ver the whole energy range from < 100 MeV to 10^ GeV. 
Our calculationsE3'c2l now include a correctireatment of the geomagnetic cutoff 
effects, as described by Lipari and StmeyEl Because of the large range of ener- 
gies they cover, these calculationsEaHH each give a consistent set of neutrino 
fluxes for simulating the full range of experimental data, from contained events 
to neutrino-induced muons at high energy. The algorithms used for these and 
other calculations of the atmospheric neutrino flux should be checked by com- 
paring their corresponding muon fluxes with the full set of measurements of 
muons at high altitude that are becoming available. 

3 Neutrino oscillation interpretation of the flavor anomaly? 

As shown in the previous section, the predicted neutrino flavor ratio is quite 
robust because many sources of uncertainty cancel in the ratio. The most in- 
tensively investigated physics explanation of the atmospheric neutrino anomaly 
is the possibility of neutrino oscillations. As an example, consider j/^ ^ i^r 
oscillations, which occur with probability 

P.„^.^ = sm'2e sin^ (^1.27 .^m^Jg^^ . (4) 

There is no visible up/down difference in contained events, for which ^ 
1 GeV. Since L ^ 20 km for downward neutrinos this gives a lower limit of 
approximately, 6m'^ > 0.005 eV^. If the effect begins to disappear for higher 
energy, there would then be an upper limit on Sm?. This is the significance 
of the Kamiokande extended analysis that includes a multi-GeV sample of 
data. In the multi-GeV sample there is an apparent up/down difference in 
the comparison between calculated and observed flavor ratio. The anomaly is 
largest for upward neutrinos with pathlength L ^ « 6000 km. The ratio 
is consistent with expectation for vertically downward events. 

The Kamiokande combined analysis defines allowed regions for both Vf, <-> 
z/^ and Vfj^ ^ Vt oscillations with large mixing angle and 5iv? in the range 0.01 
to 0.02 eV^. The preferred interpretation depends both on the normalization 
and tiie shape of the calculated neutrino flux. For the-jpiitained QCl and sub- 
GeVtSl events the calculations with high normalizationc3^c3 surest oscillations 
primarily in the v^^ <-> Vr sector whereas the calculation of BNI23 prefers <-> 

There i^an interesting energy-dependence that shows Jip in ilie multi- 
GeV eventsp^ which has been noted previously (e.g. in Refs.c3 andE3). When 
compared to the HKHM neutrino fluxes, the ratio of measured/calculated for 
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Table 2: Status of i^-induced upward muons 

IMB@ BaksanS Kamioka'^'^^Q 



Observed 



Calculated „ 
Bartol fluxL3 
Volkova fluxt 



617 



559 



1.97±0.10 



— 580 

600 — 



2.15 
2.06 



MACRO t 



255 



315 
286 



'^"^Entries for Kamiokande give fluxes in units of 10 ^"^cm ^s ^sr ^. Other 
experiments quote total number of events. 



electron-like events is ~ 1.5 as compared to ~ 1.1 for the sub-GeV sample. The 
corresponding numbers for muon-like events are ~ 0.83 and ~ 0.66. A neutrino 
spectrum that is sufficiently harder could be made to give th£_same ratios of 
measured/calculated for the sub-GeV and multi-GeV samplesEj (The statisti- 
cal significance of the anomaly would not be changed, only its interpretation.) 
Inspection of Table 1 shows that the BN spectrum has roughly the right degree 
of hardness to keep the ratio of measured/calculated constant for each neu- 
trino flavor, but at the cost of a rather extreme assumption about the nature of 
production of low-energy pions in collisions of 10 — 30 GeV protons on oxygen; 
namely, that the yield is about the same below ^ 2 GeV as for collisions on 
protons. In a recent comprehensive three-flavor treatment of the atmospheric 
neutrino anomaly, Fogli et aleS find that in fact i/^ <-> i^e oscillations provide a 
better fit to the atmospheric neutrino anomaly than t/^ <-> z^t-. 



4 Neutrino-induced upward muons 

Another way to explore higher neutrino energies is to use upward muons. Table 
2 summarizes some comparisons between experiments and calculations for the 
flux of neutrino-induced muons with energy greater than a few GeV (the exact 
value depends on the experimental cut). Tha calculationa_are shown for two 
different assumed neutrino spectra, Volkovac3 and Bartolpl and assuming no 
oscillations. The calculations also depend on the structure functions used to 
obtain the cross section for v + N ^ fi -LX. 

The 1MB calculation is an early result ll which uses the EHLQ@ structure 
functions which probably give too low a value oLlhe cross section and hence 
underestimate the expected rate by some amountEEl. BaksancZl and MACROcfl 
bothj-calculate the cross section using Morfinj-^ Tung structure function Bl- 
DIS e2I, while Kamiokande El use MRS set Ge3 The results are inconclusive 
because interpretation depends on comparison with an absolute calculaton. 
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For example, Frati eL alc3 considered an oscillation purely in the <-> ly^ 
sector. Using OwensEil structure function for the cross section, they found 1.61, 
1.97, and 2.33 x IQ-i^cm-^sr-is-i for Sm^ = 0.01 eV^ and sin2(26') = 1.0, 0.5 
and 0.0 (no oscillations) respectively. 

Comparison to the muon flux is also relevant for high energy neutrino 
fluxes, though the constraint becomes less restrictive as energy increases be- 
cause at high energy (> 100 GeV) a relatively larger fraction of neutrinos comes 
from kaon-jiiftcay as compared to the muons, which are always dominated by 
pion decayE3'E3 The muon fluxes corijespDpding to several different calculations 
of the neutrino flux at high energy rir rE3 are compared with a compilation 
of measured vertical muon fluxes from 1 to lO'^ GeV in Ref. H. All these cal- 
culations show some tendency to be higher than the mesaurements from 10 
to 100 GeV and somewhat below the data fmm 100 to 1000 GeV. This small 
systematic effect, which is also noted in Ref.Ej is not presently understood. 

The angular dependence of the upward (neutrino-induced) muon flux also 
contains in principle information relevant to an oscillation interpretation be- 
cause the pathlength changes with angle. There are, however, significant sys- 
tematic uncertainties because the experiments in general have acceptances that 
depend on direction. |(See for example the discussion of the angular-dependence 
of the MACRO datalla) There is also some difference in angular dependence 
among the calculations. In this situation it might be useful to carry out a 
two-dimensional analysis of the compaxison between data and calculation, fol- 
lowing the example of Fogli and Lisi for contained events. The variables 
might be total rate and the ratio of "horizontal" /"upward." 

Concerning angular dependence, it is interesting to note that the recent 
Kamiokanda data set (364 events) □ does not fit any calculation as well as the 
earlier setc3 with poorer statistics (252 events). For example, the i'^ <-> i^r 
oscillation referred to above with Sm^ — 0.01 gives reduced values of 2.9, 
2.6, and 3.2 respectively for sin^ 26* — 1.0, 0.5, and (no oscillation). The 
corresponding values for the earlier data set were 1.9, 1.1 and 2.0 (no 
oscillation) . 

Another way to remove some of the model-dependence from a comparison 
between expectation and observation of neutripD-induced muons is to compare 
the ratio of stopping to throughgoing muonsH The stopping muons depend 
significantly on the neutrino cross section in the few GfV region, which is 
below the deep-inelastic scattering regime. Lipari et alLj have pointed out 
that the cross section is poorly known in the low energy region. They make a 
careful evaluation of the cross section in the resonance region and find-that the 
exclusion region is less restrictive than originally estimated in Ref. e3. There 
is almost no overlap between the excluded region they find and the "allowed" 
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region for oscillations in the <-> i^r sector found by Kamiokandel 



5 Conclusion 

Results from the two large water detectorsil are consistent with each other, 
and they report a significant anomaly in the ratio of e-like to /i-like events 
as compared to what is expected if all the events are due to interaddoiis of 
atmospheric neutrinos inside the detectors. Iron (tracking) detectors t2nld are 
consistent with no anomaly but with lo}sz^|Slatistics. Preliminary results of 
the Soudan experiment are intermediatcE3ilj Calculated ratios of z^e/^'/i are 
the same within ±5% in all calculations. Main sources of difference in overall 
normalization and shape amoag the calculations of the flux oLMmiispheric 
neutrinos have been identificdEJ There are updated calculations l3CjO which 
cover the whole neutrino energy range from < 100 MeV to 10** GeV_ 
work is needed to compare these calculations with measurements I 
muons in the atmosphere at various cutoffs, times and altitudes. 

The angular dependence of the Kamiokande multi-GeV eventsPJ in com- 
bination with their sub-GeV eventsP suggests a limited range of oscillation 
parameters with large mixing angles and Sm^ in the range ^ lO"'^ eV^. Be- 
cause they require comparison to an absolute calculation (rather than a ratio 
as for lJae_mnljained events), the measurements of neutrino- induced upward 
muons cZrEaES'Ll are inconclusive at present. 

First results from Super- Kamiokande are expected very soonE3 Because 
of its large fiducial volume Super-K will be able to accumulate of order lO'* 
events in a few years. This is enough to give ~ 100 events from each cone 
of half- angle 10° in the sky. Super-K should therefore be able to pick out 
directions with particularly high and particularly low values of geomagnetic 
cutoff and therefore demonstrate that they can see the appropriate changes 
of rate characteristic of neutrinos produced by cosmic rays in the atmosphere. 
Similar remarks can be made about effects of solar modulation as we go from 
the current epoch of minimum solar activity into the next solar maximum in 
1999. In addition the systematics will improve as higher energy events will 
be fully contained. 
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